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Abstract: We theoretically study the deterministic generation of photon 
Fock states on-demand using a protocol based on a Jaynes Cummings 
quantum random walk which includes damping. We then show how each 
of the steps of this protocol can be implemented in a low temperature 
solid-state quantum system with a Nitrogen- Vacancy centre in a nano- 
diamond coupled to a nearby high-Q optical cavity. By controlling the 
coupling duration between the NV and the cavity via the application of 
a time dependent Stark shift, and by increasing the decay rate of the NV 
via stimulated emission depletion (STED) a Fock state with high photon 
number can be generated on-demand. Our setup can be integrated on a chip 
and can be accurately controlled. 
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1. Introduction 

The generation of on-demand photonic Fock states is at the heart of many photonic quantum 
technologies. Single-photon sources have been realized in a variety of quantum systems such as 
nitrogen- vacancy (N V) centres in diamond IjJ , or using quantum dots ^ . However the creation 
of photonic Fock states with a high photon number is an open challenge to date. In this paper 
we propose a novel Jaynes Cummings quantum random walk protocol that drives the cavity to 
accumulate a photonic Fock state deterministically. We describe in detail how to implement the 
theoretical protocol using a Nitrogen- Vacancy defect in a nano-diamond evanescently coupled 
to circulating light modes in a high-Q toroidal resonator at moderately low temperatures (< 10 
°K). We show that even in the case where one has error in the timing of the control pulses there 
is very high probability to generate photonic Fock states up to n = 6. 

Synthesising high-number Fock states has received much attention in the literature. A high- 
number Fock state has been conditionally produced with a probability P„ via the state collapse 
from a coherent or thermal state |l3]|4]|5l. The probability P„ of success is equal to the initial 
overlap probability of the target Fock state with the initial state of light (P(n) = Tr[|«)(«|p,„„]), 
and this probability can be quite low: for \n = 3), ~ 0.22 if starting from a pure coherent 
state \a - V3). An arbitrary quantum state of a cavity field can also be engineered if the cavity- 
qubit coupling can be very accurately controlled and recent experiments using low-temperature 
superconducting circuit-QED have synthesised microwave cavity states up to nine photons lH 
HIHl. When an excited two level system interacts with a cavity mode via the Jaynes-Cummings 
(JC) interaction of strength g, the emission probability Pemiiig, t, n) of the excited atom depends 
on the duration of coupling r and the choice of Fock state n. For certain values, terming trapping 
values, of {g,T,n), Pemit vanishes and a Fock state can be trapped in the cavity ||9] [lOl [TT| . In 
this way by sending a train of Rydberg atoms through a superconducting high-Q microwave 
cavity, Walther et al. trapped a microwave Fock state ifTOl [TTl . Indeed all the experimental 
demonstrations in high-Fock number state generation have been in the microwave regime (S] 
[TOlflTI . At optical frequencies. Brown et al. propose a system of three-level atoms in a high- 
finesse cavity ||T2) . for Fock state generation but this requires the preparation of complicated 
nonclassical states of the atoms. Until now only single photon sources at optical frequencies 
have been realized in solid-state quantum systems. 

Jaynes Cummings Damped Quantum Random Walk:- A coined quantum random walk 
(QRW) involves a coin, which we take as a qubit with Hilbert space 'He - span{\e} ,\g}}, to- 
gether with a walk on the discretised non-negative real line 1-1^ = span{\n};n = 0, 1, ■ ■ ■ ). The 



normal coined quantum random walk on the full real line (-00 < n < 00), is an iteration of a 
basic step involving a conditional displacement of the walker on the line depending on the in- 
ternal state of the walker Ud = |e){e| ® |« + l)(n| + \g){g\®\n- l)(n|, followed by a "scrambling" 
of the internal state of the walker by the action of a Hadamard operation on the internal states. 
This coined version of the QRW where the walker moves on the discretized real line Z has been 
studied intensively over the past decade. In the following we will examine the case when the 
space upon which the walker walks is the Fock ladder, « e Z*, i.e. the non-negative integers. It is 
no longer possible to have a unitary operator that implements a conditional displacement with 
constant displacement independent of the position of the walker To achieve a unitary operation 
for the conditional displacement the walker can execute a step up/down the half-line with "step 
sizes" that depend on n. Our QRW step will consist of a period of Jaynes Cummings (JC) evo- 
lution between the internal states of the coin and conditional displacements up/down the Fock 
ladder, followed by a manipulation of the internal states of the walker. Rather than a complete 
scrambling of the internal states we will just consider a flip where \g} <-» |e), are swapped. We 
have found that such a unitary QRW on the half line using the JC walk step exhibits complex 
temporal dynamics but simplifies greatly when we allow for periodic damping of the internal 
state of the coin. 

We now consider analytically the above QRW on the half-line and derive a formulae for the 
resulting map on the reduced Fock space of the walker. We will see that, if starting at \n - 0), 
the walker will, on average, step to greater values of n and will hit a ceiling value of n that 
depends on the chosen value for the JC interaction strength/time. Let the JC Hamiltonian in the 
RWA be given by Hjc - g(\e) {g\^a + \g) (e| then the resulting unitary evolution operator 
Ujc{t) - e~'^-"^'^, can be expressed as 



where N-a'a is the photon number operator. Considering the initial product state for the 
density matrix of the coin and the walker to be pc ®pw, then following evolution by the JC 
Hamiltonian we obtain pc^Pw Ujcipc® Pw)U^jq- Subsequently we allow spontaneous 
emission (amplitude damping channel) to operate on the atomic system i.e. 



where idw stands for the identity map in walker's space. Here &se — AdSQ+ AdS\ is the 
spontaneous emission channel with non-unitary Kraus generators 



where 77(f) = e~'' a positive parameter quantifying the degree by which the atomic system is 
reset by the channel, with t the nominal time over which the channel operates and T a constant 
characterising how rapid the reset process it. We have also used the notation of the adjoint 
action Ad{A) of an operator A on some other operator X as follows: X — > Ad{A)X - AXA\ 
noticing the property Ad{AB)X - Ad(A)Ad{B)X. For a general pure input state of the coin; 
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In view of the limit lim,^oo rjit) - 0, and normalization relation lap + = 1 , the last expression 
leads to the reset state \im,^oo &s E{pat](t) = IgXgl - 

Next we form the composite map beginning with the JC unitary, spontaneous emission chan- 
nel and then finally a flipping of the atomic state X = exp(-in/2o-x) and denote the entire process 
by 6: 

&= AdXo{£,sE^idw)oAdUjc, (6) 

where S acts in total coin-waUcer (atom-mode) density matrix. Choosing the action £(|e)(e| ® 
pw), gives in], 

£(|e><e|®pw) 

' Swipw)- COS(gT yfW+l)pw COS(gT VA^+1)77 /fl^ sin(,^rg4) ^^ ^^^^^^ yjN+l) y/lj " 
-icos{gT^/NTT)pw'-^^^^^^^a^fii cosigT ^/NTT)pw cosigT y/N + T)Tj j' 

(7) 

where the positive map 

; ; +sin(pT VA^h- 1) sin(eT VJvTT) 

Sw(pw) = cos(gT VA^ + l)pw cosigT ylN+l) + a^ — ^ , pw — , (8) 

yN+l yN+l 

appearing above can be shown to be trace preserving i.e. 

Tr£w(Pw) = Tr{[cos(gT VA^-h1)' + sm(gT yfN + Tf]pw] = Trpw = 1- (9) 

The model described leads to a sequence of walker's density matrices flA], and after m > 1 
steps the reduced state of the walker is 

p<;'> =Tr,[fi'"(|e><e|®pw)] =£w(pw)+a'7^). (10) 
Consider the case of number state input for the walker pw - \n} {n\. Then 

Sw{\n} <n|) = cos(gT ^n+lf |«> <«| + sin(gT V« + 1)^ |n + 1 )<n + 1 1 , (11) 

and repeated action of S in the rj ^ Q limit leads to a progressive increase in Fock number «. 
From the form of Sw in this limit we now observe that we can halt this upwards motion to 
accumulate at the trapping value n = nj, if we choose the Jaynes Cummings coupling strength 
and duration t to satisfy gr y/iir + 1 = kn, where keZ. This is the main result of this section: by 
executing a sequence of operations: Jaynes Cummings for a period, followed by spontaneous 
decay and then a complete flip of the coin space, one can, with unit probability, arrange for the 
walker to reach a steady state at the position n-nj- The behaviour of the quantum random walk 
on the half line with and without damping can be observed in Fig.[T] The dramatically different 
behaviour from the position independent step operator used in the conventional quantum walk is 
apparent. Next we show how to implement this map for an optical cavity-QED setup consisting 
of a Nitrogen- Vacancy centre in a nano diamond coupled to a high-Q optical cavity to produce 
multi-photon optical Fock states. We perform more detailed numerical simulations taking into 
account cavity and atomic decay to determine a figure of merit to synthesis Fock states of hght. 

Implementation:- To implement the above protocol we propose to use a single nitrogen- 
vacancy (NV~) center in a nanodiamond coupled to a high-finesse toroidal optical cavity at 
low temperature, while the latter is also connected to an optical interferometer and where the 
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Fig. 1. (Color Online) Jaynes-Cummings quantum random walk: Plots showing how the 
walker evolves starting in the vacuum, i.e. Trj|;3>{/?|£™[|e>{e|0|O>(O|]), as a function of the 
number of steps m and Fock number n. (a) for a completely unitary evolution with fi = 
AdHoAdU jc, where one executes a Hadamard on the coin space (b) completely unitary 
evolution with fi = AdXoAdU jc, where one executes a n flip instead of the Hadamard and 
(c) including spontaneous damping channel £ = AdXo Ad&sE °AdU jc, and the Jaynes 
Gumming coupling strength and duration chosen so that \n = 16) is a trapping state. We see 
that the latter evolution clearly leads to accumulation of the walker at the target trapping 
state. 

NV's optical transition is initialised via optical pumping, brought in/out of resonance with the 
cavity via Stark shift tuning resulting from an electric field, and undergoes periodic optical n 
flips via resonant optical laser pulses. In more detail: when the cavity interacts on-resonance 
with the zero-phonon line (ZPL) of the single NV center, the de-excitation probability of the NV 
(treated as a two level system (TLS)), [and consequently excitation probability of the cavity], 
is given by Pemit{g,T,n) - sin^( Vn + Igr), where n is the number of photons in the cavity, g 
is the JC coupling strength, and t is the interaction time. Choosing t - tt such that the nr 
photon is trapped in the cavity we have Pemit(g,TT,nT) - 0. Using a fixed tj as the time step in 
the damped JCQRW above and starting the cavity in the vacuum state leads to the cavity field 
undergoing a deterministic ratchet-like increase in Fock number until it accumulates at n = nj-. 
The trapping-state condition means that the field in the cavity reaches an upper bound and is 
prevented from being excited to a higher photonic number state. Thus via a precise control of 
the Jaynes Cummings coupling tj, an on-demand Fock state can be deterministically trapped in 
the cavity starting from the vacuum state. To do this we start the following process (Eq.|6]l from 
the excited state of the NV center, which is resonantly prepared by a .tt laser pulse: (i) we first 
switch on the JC coupling by tuning off the electric field. During this stage, the NV center emits 
a photon with the probability Pg into the cavity, (ii) After a time tt, the JC coupling is turned 
off' by bringing the NV's optical transition out of resonance with the cavity via electrical Stark 
control ifTsl and the NV center is allowed to completely decay to its ground state (GS). (iii) 
Then we resonantly pumping the NV center to its excited state again via a n pulse. Repeating 
these operations, the field in the cavity can be trapped in a selected target Fock state. 

The System:- Our setup for creation of photonic Fock state is shown in Fig |2] An optical 
toroidal cavity with high quality factor Q and resonance frequency co^ couples to the optical 
ZPL transition in an NV~ center in a type Ila nanodiamond with C^y symmetry, which is po- 
sitioned on or nearby the toroid so that it has a large evanescent overlap with the whispering 
gallery optical modes of the toroidal resonator The nanodiamond is oriented so that the [111] 
axis of NV center is fabricated to be along the z direction (see Fig O. This setup can be real- 
ized using current technology lfT6l [iTl [TSl . The toroidal cavity supports two degenerate modes, 
clockwise (CW) mode b and counterclockwise (CCW) mode a, which propagate around the 




Fig. 2. (Color Online) Solid state setup for deterministic generation of an on-demand Fock 
state of photons. A two level system (nitrogen vacancy in a nanodiamond - here shown as 
yellow at top of the toroidal resonator), interacts with counter propagating optical modes 
a and S in a high-Q toroidal resonator with intrinsic decay rate jc and which is coupled to 
a nearby waveguide interferometer at a coupling rate K^xt- Shown are the input and output 
modes a/bi„/„„i and the resulting anti/symmetric modes A„/s modes from the interferome- 
ter with each associated photon detector Dg/s- Also shown is the incident (red arrow), laser 
pulse resonant on the NV zero phonon line required to implement an optical ;r-pulse and 
the Stark shift electrode used to bring the NV's optical transition in/out of resonance with 
the cavity. Not shown is the initialising green laser and stimulated depletion laser. 

cavity along two opposite directions and form a standing wave if both modes were excited |fT9l . 
These two modes can be viewed in a basis of anti/symmetric modes A„,j - (a + h)/ V2, which 
are mutually orthogonal in space. They split in frequency due to the scattering /; from the NV 
center and rough surface. Depending on it's position the coupling of NV~ centre to one of these 
two normal cavity modes A j or A„ may occur predominantly or even exclusively with respect 
to the other mode |fT9l . In our setup, we take the NV~'s position to be at the antinode of mode 
As (node of the antisymmetric mode A„). Thus the NV center couples dominantly to the mode 
As- We neglect the small coupling to mode A^,. The mode A j is set to be red detuned to the ZPL 
transition of the NV center in the absence of any applied Stark shift ifTSl . The latter we propose 
can be used to control the coupling between the NV center and the cavity. 

An optical waveguide precisely positioned close to the cavity couples Ught in/out to/from the 
cavity with external coupling rate Kext via the input-output relations 12011211 

aout = -Oin + ^2Kext a, (12) 

bout = -bin + y/2Kexi b, (13) 

where the input and output fields of the waveguide are denoted by {ai„,bi„,a„i,t,bout}, respec- 
tively. [aM,al(t)] = [ao„t(t),aljt)] = 6(t-t') and similarly [&,„(f),^L(0] = [^««r(0,^Lf(0] = 
5{t — t'). The output fields tout and bout are mixed by a 50 ; 50 directional coupler l!22]| . We take 
the inputs to the cavity to be vacuum states, i.e. (a,„) = {bin) - 0, and thus the outputs and 
bout are proportional to a and b, respectively. Thus the outputs of the directional coupler yields 



modes Aj and A„ 1231 . leading to the detectors. Here we aim to create Fock state of the sym- 
metric mode As. Assuming the intrinsic loss of the cavity is denoted by the decay rate jc, if we 
take into account the scattering h between two modes a and b, the critical coupling condition is 
given by = ^li^+yl El. 

To switch the interaction between the NV center and the cavity, an electric field perpendicular 
to the axis of NV center is applied to induce a Stark shift. This static electric field (SEF) can be 
created by two electrodes positioned 10 jim above the setup ll24l . This distance is much larger 
than the wavelength of the field in the cavity and the extent of the evanescent field and thus 
results in negligible scattering loss of the cavity modes. During the excitation of cavity, this 
SEF is applied to shift the NV center in/out of resonance with the cavity, thus executing the JC 
step in the map Eq. (|6ll. 

A critical step in the process Eq. ^ is the rapid decay of the two level system &se{p), the 
spontaneous emission decay of the two level system (coin). This decay must be executed with a 
rate much higher than the cavity decay rate. The natural excited state lifetime of the NV ZPL is 
~ 1 Ins and this is too long to permit many repetitions of our process Eq. (|6]l even with high-Q 
cavities. To shorten this, after the JC coupling is switched off by applying the SEF we use a 
stimulated emission depletion (STED) laser beam (As ted - 775 nm), to dynamically create a 
fast decay channel from the excited state to the ground state of the NV center [|25l, and this can 
increase the effective decay rate of the NV by almost four orders of magnitude. After almost all 
of the population has decayed to the ground state \g), another laser beam on-resonant with the 
ZPL repumps the NV center from \g) to \e) ll26l[T5ll27ll28l . i.e. an optical ;7r-pulse. This is the 
final X portion of the map Eq. 

2. Detailed Model of Experimental Protocol 

We propose to implement the JC Damped QRW Fock state synthesis using a NV cavity-QED 
setup. The relevant energy level scheme for the NV center is shown in Fig. [3] (a). We neglect 
both the hyperfine electron-nuclear spin coupling and the weak electronic spin-spin interaction 
Il33l . The center has an optically allowed transition between an orbital ground state and an 
orbital excited state ^E. Both the ground and excited states are 5 = 1 spin triplets. The ground 
state has symmetry and is split into an S _^,Sy doublet 2.87GHz above an S ^ singlet due 
to the zero-field splitting 13211 . The lifetime of excited state ^E is about 1 1.6 ns, corresponding 
to a decay rate -y - 14 MHz l34l . Low temperature transform-limited single photon emission 
spectra from individual NV defects l34ll indicates that dephasing is negligible at these low tem- 
peratures. Further low temperature experiments demonstrate that the optical excited states of 
the NV can be isolated from the effects of the nearby phonon sidebands l35l at low temper- 
atures. Throughout our operation below the non-radiative decay to the intersystem state from 
\Ex,S ^) is taken to be negligible. Such intersystem crossing decay contributes to an effective 
decay to the singlet ground state 5 -. 

At low temperature (T<10 °K), strain in the nanodiamond causes the excited state ■'E to split 
into an orbital upper branch Ev and an orbital lower branch Ey (see Fig.[3](b)). Each branch is 
a spin triplet formed by three spin states S x,Sy and S -. The sublevel |E;f,5-) is well isolated 
from the other five sublevels by several GHz. Actually the state |Ey,5-) can be isolated from 
lEr,^-) because these two sublevels are associated to orthogonal transition dipoles l26l l28l . 
Therefore the spin-conserving transition 13^,5-) <-> lEy,^-), can be excited resonantly at low 
strain l27l l26l l28l . To suppress further any small spin mixing and phonon-induced transitions 
within these two excited states l26ll28l . our setup works with low-strain NV~ centers at cryo- 
genic temperatures. Moreover, the JC coupling g is assumed to be much larger than the decay 
rate y and the thermal orbital coupling and relaxation rates. In our protocol the duration when 
the NV center is excited into state |e) state is small and thus the spin mixing can be neglected 




Fig. 3. (Color Online) (a) Level diagram of a NV center showing spin-triplet ground and 
excited states, as well as the singlet system involved in intersystem crossing 1291 1301 |3T| . 
Another triplet excited state E,. is not shown here. Also shown is the JC coupling (blue) g, 
decay rate from |e) to \g), the STED illumination (red), and the laser transition for the n flip 
generated by Q ^ (orange), (b) Eigenvalues of the excited state triplet as a function of applied 
SEE 1321 . The vertical dashed line at V = marks the splitting due to the strain. At this 
position the NV" center can be excited resonantly by a ;r laser pulse flf . An electric field is 
applied to bring the NV center into resonance with the symmetric mode (resonant frequency 
(its)- (c ) Time sequence for generation of photonic Pock state showing the initialisation, JC 
coupling. Decay and X flip. 



Il26ll28l . Our protocol primarily involves the transition between |^) «-» \e), (PAt,^ -) «-> |E,,5 -)) 
and excludes other sublevels ll27l . Thus we are able to treat the NV as a two level system with 
a transition in the optical ~ 637nm. 

The dynamics of our system is given by the Hamiltonian H in the rotating wave approxima- 
tion (RWA) 

H —Hs + Hq + Hx , 

//, = - Ti(Ag + A,{t))A\As + fig[A\cr- + H.c] , (14) 
H^=m,{t)[&-+H.c.], 

where &- - \g){e\. Ag - cj^pi - cJs is the detuning of the mode A j and the ZPL transition between 
states \g) and |e) (with frequency 0J,pi), in the absence of any Stark shift. - coc + h is the 
resonant frequency of mode As shifted by the scattering /i. g is the JC coupling strength between 
a single NV center and a single photon in the cavity. In our scheme, the NV center only couples 
to the symmetric mode A.,. This is reasonable because this coupling can be predominately 
to mode As by specialized positioning of the nanodiamond jl9!|, and the scattering can also 
introduce a large detuning between the unwanted mode A„ and the ZPL of N V center The Stark 
shift As(f) is used to dynamically control the creation of cavity photons by the NV center due 
to the JC coupling ifTSl . For A 5(f) - 0, the cavity decouples from the NV center because of the 
large detuning. This process can be considered as "the JC coupling off". The excitation of cavity 
is turned on if "the JC coupling on", i.e. Av(f) = -A;,. According to our numerical simulation, a 
fast relaxation from \e) to the ground state |^) is required following the JC coupling phase for the 
preparation of a high-number Fock state with a high fidelity. Here we make use of the concept 
of "Stimulated Emission Depletion" (STED) to dynamically enhance the relaxation process 
II25I during the "decay phase" of the map (see Eq. (|6]l and Fig.[3jc )). When the STED beam is 
applied, the stimulated emission rate jsted becomes to IsTEDjIh, with Is ted /Is denoting the 
ratio of the STED pulse intensity Is ted and the saturation intensity /j. For a lifetime 1 1.6 ns. Is 
is ~ 1.85 MW cm~^ ll25l if a continuous wave (cw) STED beam is applied. A cw STED beam 
of 20 GW cm~^ can enhance the decay rate by four orders of magnitude. During the initial "JC 
coupling on" phase (see Fig.^c )), we turn off the STED beam and the nominal decay rate of 
NV center remains y ~ 14 MHz. 

We now describe the detailed steps in synthesising the process described in Eq. (|6]l. We 
assume that the optical modes in both the cavity and waveguide are initially in the vacuum state 
yielding initially zero photon number in the cavity. The time sequence for creating a photonic 
Fock state is shown in Fig[3](c): Initially the NV center is in its ground state |^) after a short 
532 nm laser pulse optically prepares the defect into the = 0(pA2,5' -)) state. A n laser pulse 
Qx is used to resonantly pump it to the excited state |e). Then the ZPL is tuned on-resonance 
with the mode As to enable the JC coupling g using a SEE. After time T7, we turn off the JC 
coupling but use the STED laser beam to create a fast decay channel to the ground state \g). 
Waiting for time Ty, almost all population decays to the ground state \g} from |e). Then a further 
optical 7T pulse generated by is applied to resonantly excite the NV center to \e) again. We 
repeat these operations until the target state is trapped. 

When Ag + As = 0, the NV center resonantly couples to the cavity mode A,. The dynamics 
of the system can then be described by a unitary time evolution operator and we further now 
assume that the time duration of this unitary may not be precisely controlled, i.e. we assume 
some noise in the target JC coupling time T7. More precisely we take Ujc - e""'?^*'"'"'*'')^'/^', 
where 6t is a normally distributed additional noise in timing with a standard deviation given by 
a parameter cr„. The excitation probability of the cavity to state |« + 1) when the NV is in the 
excited state is thus now given by Pg - sin^lgTT Vn +1(1 + St)]. Once the SEE is turned off, 
the STED laser pulse is switched on. The NV center is decoupled from the cavity and relaxes 



to its GS S ^. The timing noise during the decay is not considered in Ty because this damping 
process is insensitive to the timing error. The population in the excited state |3) decays at an 
effective rate jsted to the ground state |^). Such process can be described by a supperoperator 
e as ll36l 

&[p\^e\ I p, (15) 

where the superoperators are defined as Hsp - [Hs,p\, Lp - ycl2{2apa^ - ap - pa' a) + 
JstedI'2{1&-P&+ -&^p-p&^ with (T, - \e){e\ - and the density matrix p. After a time 
Ty = 5y, the GS S - is polarized more than 99% again. The flip n laser pulse generated by fiv 
(/tj: a! 637 nm) turns on successively to flip the NV center to the ES S -. We define a flip operator 
X - e-"'(i+<5^)o"v/2 ^j^jj = (T+ +cr_ to model this flip process as XpX\ 6x is a noise having 
the same statistic property with but independent of St- Then the density matrix after I + 1 steps 
is determined by a recurrence relation 

Pm^X6se[UjcpiU]c]xK (16) 

The system is initialized in the state po = |e)(e| ^p^ac where p^.^^,^. is the density matrix of vacuum 
state of cavity mode. 



3. Results 

Next we discuss the generation of the photonic Fock state. Throughout our simulation below, 
we neglect the excitation of A^,. This requisite can be satisfied by positioning the NV center at 
the antinode of A , [ 19 | or introducing a large scattering. 

To control the JC coupling, we consider a setup shown in Fig. |2] in which the cavity is 
designed to be off-resonance with the transition \g} <-> |e) such that |Ag| » A detuning of 
Ag = IQg is large enough to decouple the cavity from the NV center. To switch on the JC 
coupling, the transition \g} ^ \e} is tuned to be on resonance with the cavity, i.e. Ag + A^ = 0, 
by the Stark shift A, induced by the SEF ifTSl . However during the JC step the NV center will 
decay to its ground state \g} at a rate given by y and this decay during the JC step decreases the 
fidelity of the target Fock state. We assume a static, larger JC coupling g - BOy to suppress this 
detrimental process. By assuming a good optical cavity with y, in combination with the 
large JC coupling strength further improves the ultimate fidelity of the trapped photon state. 

In the ideal case of no cavity decay, a complete switching on/off JC coupling and no timing 
error, one can stably trap a Fock state with unit fidelity, see the black solid line (i) in Fig.|4] 
(a) for instance. 

Now we discuss the influences of timing error and the decay of cavity on the generated 
trapped photonic Fock state taking as an example nj - 6. Before the fifth step, only the Fock 
states with \n < 6) are excited, because in each step the prepared photon number state can only 
excite the next one. As the operation continues, the target state Inj = 6) is essentially populated. 
It can be clearly seen from Fig.|4](a) that the fidelity F - Trlppj] with Inj- - 6) Il37l increases 
quickly at first because the excitation probability P{g,TT,nT = 6) is large when the population 
in \n = 6) is small. As more population transfers to \n = 6), the probability to excite this state 
decreases. Overall we have found the generation of Focks states is fairly robust with the number 
of repetitions of the process Eq. (|6]l, providing the timing error is not too large (ii and iii). The 
cavity mode becomes stable after ~ 73 steps. 

If only the cavity decay is included (red dashed-dotted line (ii) in Fig.|4]i(a), the target state is 
stable once it is prepared after about 100 steps, and then the fidelity F is very high, about 0.97. 
The loss of cavity photon cancel the small probability of pumping from \n - 5) to \n - 6} if F 
is large, and consequently leads to the reduce of fidehty. Thus the excitation of state \n - 5) is 
considerable, see red bar in Fig.|4](b). 
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Fig. 4. (Color Online) (a) Time evolution of fidelities of target Fock state \n = 6). (i) 
cr„ = 0,7r = 0; (ii) cr„ = 0,7^ = O.I7; (iii) cr„ = l,y,. = 0.17%; (iv) cr„ = 2%,jc = 0. (b) Prob- 
abilities of photon number states at step 73. Red bar for only cavity decay cr,, = 0,yc = 0. 17; 
blue bar for cr„ = \%,yc = O.I7. Other parameters are g = 30-y,Ag = SOOy. 



(b) 



We notice that the timing error in the JC coupling causes a leakage of the population to 
higher photon number states. This leakage results in a reduction of the fidelity as the operation 
continues (iii and iv). To provide a limit for the fidelity of a Fock state we can prepare with nj- < 
6, we perform the simulation including both kinds of imperfection: the timing error (cr„ = 1%) 
and the cavity decay jc - O.I7. In this case, the probability of |n = 6) is about 0.9 from step 64 
to 94 (blue dashed line iii in Fig. |4](a)). Obviously, the prepared Fock state is stable within a 
wide operation step range. This is an advantage of the trapping state ifTOl . However about 6% 
population leaks to the higher photon number state (blue bar in Fig.|4](b)). 

A timing error of 2% substantially destroys the trapping condition (iv) and causes a con- 
siderable excitation of higher state. As a result, a considerable part, about 15%, of population 
leaks to higher photon number states. The fidelity of the target state \n - 6) decreases fast from a 
maximum F = 0.81 after 60 steps. However if the operation stops at step 60, one still can obtain 
the Fock state \n = 6) with high probability. A large error in the interaction time is the crucial 
reason why a train of atoms successively entering a cavity can not trap a Fock state with high 
probability IfTTllSSll . Thanks to a solid-state setup, this timing error can be much smaller during 
our operation. In contrast, our setup can generate a higher Fock state with a higher fidelity. 

To generate a Fock state \n = 6), we need jsted > 10'*7f To perform these simulations 
we used a cavity with g ~ 3 x 10^ lESl HO] ST] |42l |43l EM, corresponding to a decay rate 
of jc ~ 27ix\A MHz. Theoretical studies predict much higher Q factors are possible but we 
choose a lower limit in our example here B3ll44ll . We use the typical value jl2n - 14 MHz for 
the decay rate of the excited state |3) of a single NV center Il34l |45| . This decay rate can be 
enhanced by four orders in magnitude if a 20 GW cm~^ cw STED beam is applied. To suppress 
the decay of population from the state |3) during the JC coupling on, we need a JC coupling 
strength g — 307 ~ 400 MHz, which can be reached in the current experiments ll46ll . This large 
value of g allows for shorter tj and on these times scales the mixing between states E, and Ey 
is negligible. For this coupling strength, a Stark shift of Aj - \Qg ~ 2;7rx4 GHz is large enough 
to switch on/off the excitation of cavity. Such Stark shift can be created using two electrodes 
separated by 10 jum and positioned 10 fim above the NV center 12411 . 



4. Conclusion 



In conclusion, we have proposed a solid-state setup consisting of a single NV~ center and 
a high-Q toroidal cavity for the generation of a multi-photon optical Fock state through the 
iteration of a damped Jaynes Cummings quantum random walk. By iterating this walk step we 
found a method to trap an on-demand photonic Fock state with high fidelity within the cavity. 
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